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KINETIC THEORY OF BIMOLECULAR CHEMICAL REACTION, 

DIFFUSIVE DRAG, AND OTHER PROCESSES 
IN A GAS MIXTURE 

By Richard A. Hord 
Langley Research Center 

SUMMARY 

For two gas species with a temperature difference and a bulk velocity difference 
neither of which is necessarily small in magnitude, the kinetic theory of gases has been 
used to derive explicit expressions for the collision frequency, diffusive drag force, 
molecular translational energy transfer rate, and the bimolecular chemical reaction 
frequency. The derivations, which are based upon hypothesized mutual collision diam- 
eters, activation energies, and steric factors, are of interest in connection with theoret- 
ical studies of low-pressure gas mixtures with large departures from equilibrium. A 
binary temperature concept is introduced as an aid in condensing and interpreting the 
expressions derived from the kinetic theory. The expression derived for the diffusive 
drag force is used to give a more definite form to the equations of motion of the individ- 
ual species in a mixture of several inter diffusing gases. 

INTRODUCTION 

The theory of diffusion and chemical reaction in gas mixtures must be extended 
substantially if it is to meet effectively the challenging demands of a number of recent 
problems encountered in aerospace research, for example, the theoretical problem 
associated with the release of a gas in the upper atmosphere which reacts with the 
ambient atomic oxygen to produce chemiluminescence. The problems in question are 
characterized by low gas pressures and by large departures from thermodynamic equi- 
librium which are not adequately represented by the Chapman-Enskog theory. 

The branch of kinetic theory developed largely by Chapman and Enskog (see, for 
example, refs. 1 and 2) provides a rather detailed treatment of gas mixtures which are 
locally near equilibrium, but it embraces neither diffusion nor reaction kinetics when 
diffusion speeds are not small in comparison with typical molecular speeds. Large 
departures from equilibrium are, of course, admissible in certain special theories, for 
example, the kinetic theory of shock waves, that of individual or crossed molecular 



beams, and Maxwell's theory (ref. 3) of the dynamics of gases the intermolecular forces 
of which vary inversely as the fifth power of the distance between centers. 

An approach to the approximate dynamical analysis of a mixture of interdiffusing 
gases which are subject to chemical reactions and other phenomena can be based on a 
mathematical model in which each gas species has a molecular velocity distribution 
which is locally Maxwellian with respect to a reference frame moving with the species. 
Some situations of gas mixtures in which such a distribution of molecular velocities is 
especially appropriate have been discussed by Grad (ref. 4), Morse (ref. 5), and others. 
Before the accuracy and usefulness of this kind of model in a given nonequilibrium situ- 
ation is tested theoretically or experimentally, it is of interest to ascertain various 
implications by performing appropriate kinetic theory derivations. Representative 
derivations of this type and the statement of the results in forms which should facilitate 
their application and the understanding of their physical significance comprise the main 
content of this paper. As an example of the application of these results, the expression 
derived for the drag or tractive force acting between two gas species is used to formu- 
late equations of motion of the individual species in a mixture of several interdiffusing 
gases. The dynamical model envisaged encompasses diffusion speeds which are not 
necessarily small compared with thermal molecular speeds. 

When the present paper had been completed, the author learned that some of the 
kinetic theory results had been given previously by Tanenbaum (ref. 6), who also con- 
sidered some other intermolecular potentials and compared the results with those of 
Maxwell (ref. 3). 


SYMBOLS 

b impact parameter (see, for example, ref. 13) 

Ci mean molecular velocity of ith species (ref. 1); equivalent to u^ of present 

paper 

Dj^ drag force, per unit volume, exerted by species j on species i; 

magnitude, 

sum of drag forces on species i 
E„ activation energy per reactant pair of molecules 

a. 

E a kinetic energy of molecular pair relative to their center of mass, along line 

of centers before collision 
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body force, per unit mass, on species i 
divergence of stress dyadic of species i 
I impulse on molecule 2 in collision (see eq. (C3)) 

1^ defined by equation (E6) 

Ig defined by equation (E9) 

Kjg defined in appendix C 

M(v) defined by equation (A14) 

mass of molecule of species i 
m.m. 

m.. reduced mass, - 1 — 

i] n^ + nij 

N(v) defined by equation (A3) 

Nj 2 frequency of collisions, per unit volume, per unit time, involving one 

molecule of species 1 and one of species 2 

t 

Njg frequency of collisions, per unit volume, per unit time, involving one 

molecule of species 1 and one of species 2 and leading to chemical 
reaction 

n^ number density of species i (molecules per unit volume) 

n summation index 

P(v) defined by equation (B2) 

P defined by equations (El) 

P ^2 probability or steric factor 

Q(v) defined by equation (B4) 



defined by equation (BIO) 


R 12 


S(v) 

T,Ti 

T 12 

t 


u* 


V 


V V V 
v x ?v y 5 v z 


v 


V 



dimensionless relative speed defined by equation (3) 
defined by equation (CIO) 
integration variables 

absolute temperatures; subscript i refers to species 

binary temperature defined by equation (10) 

time 

bulk velocity of species i 

defined by equations (A4) or (A5) 

rectangular Cartesian components of V 

magnitude of v 

bulk relative velocity, Uj - Ug 

molecular velocity for species i 

magnitude of v^ 

rectangular Cartesian components of v^ 


i 


W activation speed defined by equation (E2) 

W(v) defined by equation (C13) 

w magnitude of w 

w j magnitude of v^ - v 
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z 

a 


01>/3 2 >0i2 

°12 

e 

0 w 


ij 


X 


V 

p i 


magnitude of Vg + v 

molecular relative velocity, v^ - Vg 

rectangular Cartesian components of w 

molecular approach speeds defined by equations (Dl) 
dimensionless activation speed, 

dimensionless parameter defined by equation (2); see figure 1 
defined by equations (4) 

diffusive drag coefficient (ref. 1, p. 415) 
co-latitude in spherical coordinate system 
angle between a and w immediately before collision 
dimensionless parameter defined by equation (B14) 
defined by equation (B14) with r replaced by r^ 
vector defined by equation (Cl); magnitude, X 
number of gas species in mixture 
density of ith species 


a i + °\ 

ct.. collision diameter, 1 

!] 2 

diameter of molecule of species i 

energy flow from species 1 to species 2, per unit volume, per unit time 



u 12 dissipative heating rate of species 2 due to diffusive drag of species 1, per 

unit volume 


(p azimuth in spherical coordinate system 


cp w angle defined following equation (C7) 

u> dimensionless parameter defined by equation (11) 

£ dimensionless relative molecular speed, ? 12 1/2 w 

V gradient operator 


f 


p oo p a 

\ • • • ) 

_ OO ^ c 


An overbar denotes a mean value of associated quantity. An arrow denotes a 
vector. An underline denotes volume element in velocity space, for example, 


dV = dV x dV y dV z . 


ANALYSIS 

The kinetic theory of gases will be used to derive analytic expressions for those 
properties or rates which are most essential to a dynamic description of a chemically 
reacting mixture of gases with different temperatures and different bulk velocities. 

In the mathematical analyses performed herein, each gas species is assumed to 
have a velocity distribution which is locally Maxwellian when viewed from an inertial 
frame moving with the bulk velocity of the species. The smooth, elastic, hard-sphere 
molecular model is employed because of its simplicity, its use in yielding reference 
values to which more refined values can be normalized (see ref. 2, pp. 523ff), and the 
fact that it roughly approximates the strong, short-range repulsive forces which often 
predominate in molecular collisions. Since diffusion speeds are not restricted in the 
present study and temperature differences between species are admissible, this analyt- 
ical work represents an extension of the existing kinetic theory of gas mixtures of rigid- 
sphere molecules. (See refs. 1, 2, and 7.) Expressions in closed form are derived for 
the interspecies (binary) collision frequency (appendix A), diffusive drag (appendix B), 
molecular translational energy transfer rate (appendix C), definition of a binary temper- 
ature (appendix D), and number of bimolecular collisions resulting in chemical reaction 
(appendix E). In appendix E, the additional assumptions required, which are simple and 
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well known, are concerned with the nature of molecular collisions leading to chemical 
reaction, the quantum theory of which is largely undeveloped. (See, for example, ref. 8.) 

The expressions derived herein are asymptotically equivalent to well-known 
expressions of kinetic theory for equal temperatures and relative species velocities 
which are of small magnitude compared with typical molecular speeds. 

Assumptions 

The detailed kinetic -theory derivations given in appendixes A to E of this paper are 
based upon the assumptions that: 

(1) Binary molecular collisions prevail. 

(2) Molecules collide as smooth, elastic, hard spheres. 

(3) Each gas species has a Maxwellian velocity distribution when viewed from an 
inertial frame moving with the species bulk velocity. 

(4) Collisions leading to chemical reaction are delimited by a steric factor and an 
activation energy. (This assumption is needed only in appendix E.) 

Temperature and bulk velocity are both permitted to differ, without restriction, from one 
species to another. No explicit consideration is given to ionization, dielectric polariza- 
tion, or diamagnetic or paramagnetic susceptibility, even though the results may be of 
some value in studying problems involving these phenomena. 

In connection with assumption (2), it will be seen that only the mutual collision 
diameters, not the individual radii, appear in the final equations. Moreover, it will be 
observed that the notion of collision diameter is interpretatively flexible in the case of 
the chemical considerations of appendix E, because it occurs multiplied by a steric or 
probability factor. 

In the dynamical considerations of diffusion in gas mixtures, no refinement has 
been made to account for thermal diffusion. It is assumed here that the interaction 
between the gas species in a mixture is solely that of the diffusive drag forces associated 
with species bulk velocity differences. 

Principal Derived Expressions 

In this section, the principal expressions which are derived in the appendixes are 
presented. For additional discussion and details the reader is referred to the appropri- 
ate appendix. 

The expression derived in appendix A for the number Njg of collisions between 
molecules of species 1 and species 2 per unit volume, per unit time, can be written 
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( 1 ) 


N 12 = (87 r) 1 / 2 !!.,]! 


l ll 2°12 




m i 




) 


where n^, and (i = 1,2) are the number density, molecular mass, and absolute 
temperature of species i, a ^ is the collision diameter, k is Boltzmann's constant, 
v is the magnitude of the species relative bulk velocity, and a is the nearly constant 
parameter defined by 



The dimensionless speed parameter r in equation (2) is defined by 


r 


" %2 


^v 


where, by definition, 


1 


^12 



_1_ 

02 


0i 


2kT i 

nn 




( 2 ) 


(3) 


(4) 


Equations (A19), (A21), and (A28) are, in view of the definition of a, alternative forms 
of equation (1). For v = 0 (or, equivalently, r = 0) and = T 2 , equation (1) reduces 
to the well-known collision-frequency equation of the kinetic theory of gases. (See, for 
example, ref. 1, p. 90, eq. (4).) 

The diffusive drag force D^ 2 , per unit volume, exerted by species 1 on species 2, 
is derived in appendix B and can be written in the form: 


5 12 = §<8 T ) 1/ V2^ 2 (^ + ^ + f v2 ) m 12 7 


kT< kTr 


\ 1/2 


( 5 ) 


m l m 2 _ 

where m 19 is the reduced mass , v is the bulk velocity (magnitude v) of 

L,L m 1 + m 2 

species 1 relative to species 2, and k is the nearly constant parameter defined by 
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( 6 ) 



With this definition of k, equations (B9) and (B16) are equivalent forms of equation (5). 
For small v (or, more precisely, for r « 1) and = T 2 , equation (5) yields an 
expression for the resistance opposing diffusion which is asymptotically equivalent to 
e l2 n l n 2^1 ” ^2) re ^ erence 1) page 415. The same limiting result was obtained by 

Stefan (ref. 9) and, except for the factor 4/3, by Maxwell (ref. 7), who first analyzed the 
resistance opposing diffusion. For large v (that is, for r » 1), the magnitude of the 
diffusive drag force given by equation (5) is proportional to v^. By comparison, 
Maxwell's analysis (ref. 3) of a gas mixture in which the intermolecular forces are 
repulsive and vary inversely as the fifth power of the distance between centers yielded 
a diffusive drag proportional to v, even for large v; it should be recognized that the 
analysis in reference 3 holds for all relative velocities because of the mathematical 
simplification wrought by the selection of the inverse fifth power model. 

The time rate of change of the molecular translational energy, per unit volume, of 
the second gas species (viewed from an inertial frame moving with species 2) is derived 
in appendix C and can be expressed as 


4m 


K 12 = 


12 


m l + m 2 


k( T l' T 2) N 


m. 


12 


vD 


mj + m 2 


12 


( 7 ) 


where N^ 2 is the collision frequency given by equation (1) and D^ 2 is the magnitude 
of the diffusive drag force given by equation (5). (Eq. (7) is equivalent to eq. (C17).) If 
the first and second terms on the right-hand side of equation (7) are denoted by t^ 2 and 

uj2> respectively, it is seen that is the flow of energy from species 1 to species 2 

due to the temperature difference and is the dissipative heating rate of species 2 

associated with the diffusive drag. It is also noted that 


t 12 “ -t 21 


( 8 ) 


and since Dj 2 = Dg^, 


u 12 + U 21 = vD 12 


( 9 ) 


Equation (9) states that the total dissipation rate due to diffusive drag is vDj 2> 
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In equation (94) of reference 3 (which should be corrected for typographical omis- 
sions by inserting the factor 1/2 before the expression in braces on the left side and the 
factor /3p different in meaning from of the present paper, before the expression 
in brackets on the right side), the expression on the right-hand side is analogous to the 
present one for except for the number of molecular degrees of freedom involved 

and the independence of velocity distribution made possible by Maxwell's assumption of 
an inverse fifth-power intermolecular repulsive force. If Maxwellian velocity distribu- 
tions are assumed for the two gases, the terms of Maxwell’s expression which corre- 
spond to the terms r 12 and u 12 are > respectively, proportional to Tj - T 2 (with 

no dependence on temperatures in the coefficients) and proportional to v^. Thus, the 
present assumption of elastic, hard-sphere molecules leads to a more complicated 
functional dependence on Tj, Tg, and v; see equation (7) in the light of equations (1) 
and (5). 

In the collision theory of chemical reactions (see, for example, ref. 10), the most 
significant single variable seems to be the total kinetic energy along the intermolecular 
line of centers immediately before collision. Its hypothesized threshold value is then 
called the activation energy. With the added notion of a probability or steric factor, the 
probability of chemical reaction resulting from molecular collision can be expressed in 
mathematical terms. This theoretical approach to the kinetic theory of gas reactions, 
which is of course not without shortcomings, is followed in appendixes D and E. 

For the purposes of this paper, it is convenient to introduce a binary tempera- 
ture Tj 2 > which is a measure of the collisional energy available to produce bimolecular 
gas reactions, according to the definition 


kT 12 _ kT l 


m 


12 


k ^2 2 

+ — - + x 1 v 

m< m 0 3 


where w is a nearly constant parameter defined by 


w = 


U 

4 r 2 




( 10 ) 


(ID 


The analytical basis for the definition of the binary temperature T^ is provided in 
appendix D and equations (D13) and (D14) given therein are equivalent to equation (10). 

It should be noted that the right-hand side of equation (10) is similar to forms occurring 
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in equations (1) and (5). Therefore, two other binary temperatures an <3 T^ 

can be defined by the equations 


kT 

m 


( 1 ) 


12 

12 





(12) 


kT (2) 

Ki 12 

kT x 

+ ^2 

m 12 

3 

m 2 


(13) 


and used to write equations (1) and (5) in simpler forms more akin to existing kinetic 
theory; although these equations are not needed herein, they relate to a considerable 
extent the principal derived expressions. The parameters a, k, and c o, which are 
plotted in figures 1, 2, and 3, respectively, all approach unity asymptotically as v 
(or r) approaches zero. 



.1 .51 5 10 50 100 

Dimensionless speed parameter, r 

Figure 1.- Coefficient a required for equation (A21) to be exact formula for collision frequency (and, 
correspondingly, for eq. (A27) to equal each side of eq. <A20)h 
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Dimensionless speed parameter, r 

Figure 2.- Dimensionless coefficient k plotted as a function of the dimensionless speed parameter r. 



Dimensionless speed parameter, r 

Figure 3.- Coefficient o) required for equation (DIO) to be exact. 

The expression derived in appendix E for the number of bimolecular colli- 

sions per unit time, involving one molecule of each gas species and leading to chemical 
reaction, can be written 


N 12 ~ P 12 n l n 2 cr 12 


'W /2 J_ 


, 012 / 


2r ' 


z + r 


,(z-r Y 


z - r 

e ( Z «) 2 


+ 




( 14 ) 


12 



where is the probability or steric factor and 


a \ 1/2 

\ m 12 / 

in which E a denotes the activation energy; z, like r, is dimensionless. Equation (14) 
is equivalent to equation (E14). 

If the parameter z is not greater than approximately unity, table I shows that 
equation (14) can be written in the simpler, approximate form 


. . - E a/ kT 12 

M - p XT p a / 

^12 - i i2 JN 12 e 


(15) 


This form (or the equivalent, eq. (E19)) becomes an exact representation of equation (14) 
when the parameter lo in T 12 (of eq. (10)) is replaced by the parameter u>' pre- 
sented graphically in figure 4. 


TABLE I.- RATIO OF APPROXIMATE VALUE OF n’ 12 (EQ. (E19)) TO EXACT VALUE (EQ. (14)) FOR 
VARIOUS COMBINATIONS OF DIMENSIONLESS PARAMETERS r AND z 





Ratio of N ^2 

to exact value for 

r of - 





0.1 

0.2 

0.5 

1 

2 

5 

10 

20 

50 

100 

0.1 

1.0000 

1.0000 

0.99992 

0.99950 

0.99908 

0.99966 

0.99990 

0.99998 

1.0000 

1.0000 

.2 

1.0000 

.99999 

.99970 

.99757 

.99637 

.99866 

.99962 

.99990 

.99998 

1.0000 

.5 

1.0000 

.99994 

.99833 

.98903 

.97852 

.99167 

.99762 

.99938 

.99990 

.99998 

1 

1.0000 

.99987 

.99604 

.97149 

.93050 

.96769 

.99053 

.99753 

.99960 

.99989 

2 

1.0001 

1.0010 

1.0211 

1.0642 

.93334 

.88625 

.96325 

.99021 

.99841 

.99960 

5 

1.0032 

1.0623 

2.5804 

27.630 

104.10 

1.3461 

.81726 

.94198 

.99012 

.99751 

10 

1.0688 

2.0527 

1276.8 

* 

3k 

3k 

2.4238 

.81088 

.96165 

.99014 

20 

1.9167 

171.55 

* 

* 

3k 

=k 

3k 

4.7581 

.86471 

.96160 

50 

13902 

* 

* 

* 

3k 

3k 

3k 

3k 

11.915 

.80880 

100 

* 

* 

* 

* 

3k 

3k 

3k 

3k 

3k 

23.895 


*Denotes extremely large number. 
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Dimensionless speed parameter, r 

Figure 4.- Value of of plotted as a function of r for various values of z. When substituted for 0 ) in 
equation (E19), of renders the equation exact. 


APPLICATION TO DIFFUSION IN GAS MIXTURES 

The dynamics of a mixture of interdiffusing gas species with emphasis on the drag 
forces due to substantial relative bulk velocities is considered briefly. The effect of 
chemical reactions on the species bulk acceleration, which may not be negligible in the 
case of fast reactions (for example, many reactions involving free radicals), is not con- 
sidered here. 

The equation of motion of the ith species (i = 1, . . .,v) in a gas mixture can be 
written (see refs. 11 and 12) 

< 16 > 

where and u^ are the density and bulk velocity, respectively, of the ith species, 
t denotes the time, V is the gradient operator, is the body force per unit mass 

acting on species i, Gj is the divergence of the stress dyadic (in a simple case, 
is the negative of the gradient of the partial pressure p.), and Dj is the total force per 
unit volume exerted by the other gas species on species i. If the frame of reference in 
which equations (16) are employed is not an inertial frame, each body force will 
include one or more inertial forces, for example, centrifugal force, Coriolis force 
(dependent upon u-J, etc. 

For the purposes of this analysis, the density of the gas mixture is assumed to be 
sufficiently low for binary molecular collisions to prevail; each Dj is then a vector sum 


14 



of binary drag forces, * 


that is, 




( 17 ) 


where D,, is the drag or tractive force, per unit volume, which species j exerts on 
J 1 ^ 

species i; for j = i, Djj = 0, by definition. 

The calculation of in a specific case is generally difficult. One notable 

exception is the binary gas mixture considered by Maxwell (ref. 3) in which each mole- 
cule repels every other molecule with a force along the line of centers of mass and 
inversely proportional to the fifth power of the distance between these centers, the pro- 
portionality constant being determined by the particular type of molecular pair involved; 
in this case a remarkable analytical simplification appears and is proportional to 

the product of the densities and to the relative bulk velocity of the two species no matter 
what molecular velocity distributions exist in the two gas species. Linear drag forces 
are discussed generally in reference 11. 


Another notable exception is a gas mixture whose molecules behave in collisions 
as hard elastic spheres (which is the model employed in the present paper), provided that 
the distribution of molecular velocities in each species can be adequately approximated 
by a Maxwell -Boltzmann distribution with appropriate temperature. The temperatures 
of the different species need not be equal for the calculation to be completed as in appen- 
dix B. For a mixture with more than two species, equation (B16) together with the 
assumed prevalence of binary collisions yields the nonlinear drag force 


5 jl = f<2») 


1 / 2 mij n lVlj 2 


K- , 




- \ 2 


. 1/2 


20 : 


13 


5 V 3 


“0 ( “j ■ “i) 


(18) 


Here, n^ and u^ are the number density and bulk velocity, respectively, of the ith gas 
species, cr- is the collision diameter of the molecular pair (i,j), is given by equa- 

tion (B14) with r replaced by / 3 i j 1 ^ 2 | Uj - uj , and, by definition, 


*In a real gas mixture at sufficiently high pressures, the separability into binary 
drag forces may not be possible. Moreover, the binary forces of interaction of real 
gas species may not be pure drags, but may have components normal to the relative 
bulk velocity vectors. 
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= L 

m i 

^ " 2kT t 


m,, = 


m i m i 


i] mj + mj 


( 19 ) 


where m i and T- are the molecular mass and molecular translational temperature, 
respectively, of species i. 

Equations (17) and (18) provide an explicit expression for the drag or tractive 
force in equation (16). The selection of values of the collision diameters can 

be based upon theoretical interpretations of gas or molecular data. 

If considerable fractions of the molecules of one or more gas species are in 
excited electronic energy states, the corresponding collision diameters may require 
substantial correction from those applicable to molecules in the electronic ground 
states. 

Equations (16), (17), and (18) yield species equations of motion which reduce to 
equations (1) of reference 1 (p. 415) for near -equilibrium flows for which the species 
stress dyadic is the product of the species partial pressure and the unit dyadic. 


CONCLUDING REMARKS 


Some of the basic dynamic characteristics of a model of a reacting gas mixture 
which is not necessarily near equilibrium have been examined. The comparison of the 
results with the corresponding ones of more elementary kinetic theory has been facili- 
tated by the introduction of the binary temperature concept. 

Simple models or approximations which are analytically tractable are needed in 
preliminary investigations where a high degree of accuracy is not essential. Ultimately, 
theoretical structures of greater accuracy and detail, requiring electronic computers for 
their implementation, must be developed. 

The results presented in this study provide one starting point for further develop- 
ment of the kinetic theory and continuum dynamics of gas mixtures which are far from 
equilibrium. This development can evidently proceed along several lines, such as 


16 



the incorporation of more elaborate models of molecular collisions, the inclusion of radi- 
ation effects, and the admission of more general molecular distribution functions. 

Langley Research Research, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., August 19, 1966, 

188-37-01-02-23. 
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APPENDIX A 


COLLISION FREQUENCY FOR TWO GASES 
WITH DIFFERENT TEMPERATURES AND 
DIFFERENT BULK VELOCITIES 

Consider a mixture of two gas species with hard-sphere molecules of diam- 
eters (i = 1,2) and number densities n^. Their mutual collision diameter is, then, 

°\ + a 2 

ct 10 = — — . Let each gas have a Maxwellian velocity distribution with respect to an 

2 

inertial frame and let the velocity of the first frame with respect to the second be v. 
Let the gas temperatures and molecular masses be denoted by Tj and m^, respec- 
tively. The frequency of collisions N^> that is, the number of collisions between 
molecules of species 1 and species 2 per unit volume per unit time, is derived in this 
appendix. 

With the aid of the kinetic theory of gases (see, for example, ref. 13, chs. 2 to 5, 
especially pp. 21 and 22), the collision frequency can be written, the molecular veloc- 
ities v t being referred to the second frame, 

N 12 = n l n 2 (ir) (01*2 ) V2 f we ^ ^1 ^2 < A1) 

or, the molecular velocities Vj being referred to the first frame, 

N 12 - n l n 2 (t 2 ) f we ,3lVl ^ ^1 ^2 (A2) 

where 


f.r ...r 

_ QO _ C 



(i = 1,2) 


18 



APPENDIX A 


w 


= |w| = |v x - v 2 | = [ ( v lx - v 2x ) 2 + (v ly - v 2y ) 2 + (v lz - v 2z) 2 ] 


1/2 


W 1 = Vj - V 


w 2 = v 2 + v 


v. = v. 
1 1 


(i = 1,2) 


dv. = dv* dv. dv. 
l ix ly iz 


(i = 1,2) 


For given values of the parameters n^, Oj 2 , and (3^ (i = 1,2), the collision fre- 
quency N^ 2 is a function of the relative speed v = | v | of the two gas species. Let 


N(v) = § 


“01 w l ”^2 V 2 
we 1 1 ^ ^ dv 1 dv 2 


(A3) 


which is the sextuple integral appearing on the right-hand side of equation (Al). 

In order to evaluate the integral defined by equation (A3), it is convenient to change 
from the variables ^v pV 2 ^ to variables ^V,wj accordance with the linear 
transformation 


Vj - v = V + 


02 /— 


(w - v) 


v 2 = V 


0i + h 

0i /— — \ 

— (w - v) 


(A4) 


01 + 02 


This linear transformation has a constant Jacobian, the absolute value of which is unity, 
and has the inverse given by 


— 0i /_» 0 2 _► 

V = — (v 1 -v) + — v 9 

0j + /3gV / /3 1 + ^ 2 


(A5) 


w = Vj - v 2 
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The fact that the second of equations (A5) is consistent with the earlier usage of the 
symbol w in relation to Vj and v 2 is essential to the intended application of the 
transformation. The transformation (A4) maps one six -dimensional space (v,w) onto 
the other ^Vj,v 2 ^ in a one-to-one manner. 

Equations (A4) yield 

*l(vi - v) 2 + %v 2 2 = (i3 x + /3 2 )V 2 + /3 12 (w - v) 2 (A6) 


where 


%2 " 


^1^2 _ 
Z 3 ! + @2 


Before applying these results to the evaluation of the integral (A3), it is interesting 
to note that equations (A3), (A4), (A5), and (A6) can be written in the symmetric forms 


where 


N(v) = § 


we -Vl’ 2 -^2' 2 -^l( 5 -) 2 -» 2 ] dVi , dV2 , 


v,’ = V + 


P? — 


1 " h + h 


w 


V = v - 


% - 


^1 + ^2 


w 


— ^1 — ^2 
V = - — v/ + 1 


V V ey 

+ Z 3 2 Z 3 ! + % 


w = v 1 ' -v 2 ’ 


%V 2 + W 2 = (h + hY 2 + ^I2 w2 


v l' = v l 


01 - 


Z 3 ! + Z3 2 

- » Z 3 ! _ 

Vn' = Vo V 

2 2 h*h 




(A7) 


(A8) 


(A9) 


(A10) 


(All) 
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The molecular velocities v^' and Vg' are now referred to the natural inertial frame, 

/3jV 

which moves with velocity = — relative to the second inertial frame described 

% + ^2 

earlier. In the natural inertial frame, equation (A10) reduces, in the special case for 
which = Tg = T, upon multiplication of both sides by kT, to the usual kinetic energy 
transformation to a center-of-mass system; similarly, the velocity V, by the first of 
equations (A9), reduces to the velocity of the center of mass of the two molecules in the 
natural frame. Moreover, in this case of equal temperatures, the velocity of the natural 
frame can be interpreted as the velocity of the center of mass of two molecules, one of 
which has mass m^ and is at rest in the first inertial frame, and the other of mass mg 
and at rest in the second inertial frame; this fact follows from the preceding discussion 
since v^ = v and v 2 = 0 imply V = 0. 

v l" v ) =w l> the re l a ted linear transformation (A4) 
can be used to rewrite equation (A3) in the form 




where 


dV = dV x dV y dV z 

— = dw x dw y dw z 


Performing the integrations with respect to V , 


V , and V z yields 


N(v) = 



-/3i 2 (w-v) 2 


we 


dw 


[h * %) 


3/2 -^ 12 v 5 


M(v) 


where, by definition, 


M(v) = ^ 


we 


/3i 2 (w 2 -2v • w) 


dw 


(A12) 


(A13) 


(A14) 
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Changing to spherical coordinates (w ,9,cp), with the half -line whose co-latitude is 0=0 
having the direction of v gives 


n2ir p 77 p°° 2/3 19 vw cos 0 -/3iow2 9 . , 

M(v) = \ \ l we ^ e 1Z w^sin 0 dw d0 dcp 

Jo JO Jn 

f 

@12 V J 0 


2,r ^ w 2 e ^ 2 sinh^2j3 12 vwjdw 


2tt 


f £ 2 e~^ sinh(2r£)d£ 

Jo 


% r 


where ? = 31111 r = ^12 by definition. Hence, 

M(v) = ■■■■ e f2 f ^ 2 [e‘^" r ^ 2 - e "^ +r ^d? 
^12 r J 0 


77 e r 


%2 2r 


f (s + r) 2 e" s2 ds - f (s - r) 2 e" s2 ds 

_v» *^T* 


(A15) 


(A16) 


wherein the substitutions £ - r = s' and £ + r = s" were made and the primes were 
subsequently dropped. Equation (A16) can be written as 


(A17) 

The first of the four quadratures requires one integration by parts in order to reduce 
equation (A17) to the form 


M(v) = 


r 


f (s 2 + 2rs + r 2 )e" s2 ds + f 4rse" s ds 

t/r 


M(v) = 


3 12 


e r re” r + (l + 2r 2 ) ^ e -s ds 


(A18) 


which can be evaluated numerically with the aid of tables of the error function. 
Equations (Al), (A3), (A13), and (A18) yield, therefore, 
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2 ') £ 

Equation (A19) gives, then, the number of collisions per imit volume per unit 

time between molecules of two gases with Maxwellian velocity distributions and a relative 
bulk speed v = 0^ The formula is exact for hard-sphere molecules with collision 
diameter a \2’ Equation (A19), in the limit in which the temperature of one gas 
approaches zero, reduces to the expression for the collision frequency for a molecular 
beam of zero temperature impinging on a gas at equilibrium. (See ref. 14, p. 108.) 

In equation (A19), the function (of r) in brackets is an entire analytic function 
(with removable singularity) whose power-series representation is given by the right- 
hand side of the equation 



N 12 = n l n 2 CT 12 2 f^) 6 + [l 


e“ s2 ds 


3 


(A19) 


For real numbers r ^ 0, the first two terms of the series are positive and subsequent 
terms alternate in sign. 

It will also be shown in this appendix that equation (A19) can, with both justification 
and advantage, be written in the form 


N 12 - n l n 2 a 12 2 ^j 2 ( 1 + | “ r2 ) 

= 2(2-) 1/2 ni n 2 a 12 2^ + fv2y /2 


(A21) 


where a is defined by equating the right-hand side of equation (A19) and the first right- 
hand side of equation (A21); therefore, 



(A22) 


where, by the usual definition, 

erf r = — f e -s ^ds 
/ir ^0 
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For 


0 S r £ r max = 1, the first two terms of the power series (eq. (A20) yield the 




V 


15. The 


approximation 2} 1 + ^- ) with an error which is positive and less than r max 

power series (eq. (A20)), rewritten as a series of nonpositive powers of 1/r, is the 
Laurent expansion about the ideal point at infinity; thus, the function has an essential 
singularity there and no expansion in nonnegative powers of 1/r is possible. However, 
if consideration is restricted to real values of r » 1, it is easy to show with the aid of 
tables that the left-hand side of equation (A20) is approximated by ttV 2 j^ r + (Sr)' 1 ], 
and somewhat less accurately by The preceding results suggest combining the 

features of the function on the left-hand side of equation (A20) at small and large r in 
approximations such as 


! ( i+ H 

H- 2 ) 


1/2 

(A23) 

1/2 

1 

(A24) 

A 1 /* 


r 4 j 

(A25) 


and so forth. These approximations, in ratio to the given function (A20), are plotted in 
figure 5. 

The approximation (A23) yields, with equation (A19), the approximate collision 
frequency 


N 12 = n l n 2 t, 12 2 (^) 2 ( 1 + | r2 ) 


or, upon rearranging and substituting for [ 3 ^ an d r , 


N 


12 



n l n 2°12 



+ 




(A26) 


which has an obvious physical interpretation in terms of the thermal translational 
energies and the kinetic energy, per unit mass, of either gas relative to the other. 
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z R r2 ) 


e^ds 


frequency to exact value 



Dimensionless speed parameter, r 

Figure 5.- Comparison of approximations (eqs. (A23), (A24), and (A25)) for collision frequency with exact value. 
Equations <A23) and (A24) are of the form of equation (A27) with a = Constant. 


The graphical representations of the approximations (A23) and (A24) in figure 5 
suggest expressing the function on the left-hand side of equation (A20) in the simple and 
exact form 


1/2 


2 ^1 + | o-r 2 j 

which corresponds to the collision frequency (exact) 


N,, = 4(|) 1/2 Wl2 2 (3 f vA 

12 \3 ) 1 2 12 1 2 m t 2 m 2 2 J 


1/2 


(A27) 


(A28) 


where a is nearly constant. Equation (A28) is equivalent to equation (A21) and a, 
which is given explicitly by equation (A22), is plotted as a function of r in figure 1. The 
dimensionless parameter a, as a function of r, is asymptotically equivalent to the con- 
stant 1 for r « 1, and to the constant 3jr/8 for r » 1. 


25 



APPENDIX A 


Another way of deriving the form of equation (A21), which is less intuitive and 
synthetic but more satisfactory from an analytical standpoint, consists of squaring both 
sides of equation (A19), expanding the resulting right-hand side to the second power 
in r by Taylor’s theorem with a remainder, and demonstrating (with a plot such as 
fig. 1) the limited variability of the (second derivative) coefficient of the final or remain- 
der ^r2/2lj term. Since the expanded function is known, the final coefficient can 

be evaluated for each value of r. The linear term in r is absent from the aforemen- 
tioned quadratic expression ^for Nj^) and the latter therefore leads to an expression 
for Nj 2 of the form of equation (A21). This alternative derivation provides a firm 
theoretical basis for adopting the expression for the collision frequency in the form of 
equation (A21). 
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DRAG OR TRACTIVE FORCES FOR TWO GASES 
WITH DIFFERENT TEMPERATURES AND 
DIFFERENT BULK VELOCITIES 

The drag or tractive force per unit volume exerted by the first gas on the second 
is derived in this appendix. The kinetic theory of gases (see, for example, ref. 13, 
chs. 2 to 5, especially eqs. (2.27), (2.33), and expression (5.4)) yields, for the fore- 
mentioned drag force, 



where b and m^g denote the impact parameter and reduced mass, respectively, and 
where, by definition, 

P - - -^ w* 2 -^,, 2 

P(v) = + v- wwe ^ dv^ dvg (B2) 

With the change from the variables v^Vg to the variables V,w defined by equa- 
tion (A4), equation (B2) becomes 

P(v) . ( v • w dw = /— I_\ 3/ Y' 5l2v2 Q(v) (BS> 

j — + M 

where, by definition, 

f -/3i 0 (w 2 -2v • w) 

v • w we dw (B4) 

Changing to spherical coordinates w ,6,(p with the half -line whose co-latitude is 0 = 0 
having the direction of v gives 
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n2lt nTt n°° 9 2ft 9 VW COS 9 -ftijW^ 9 

- ' * ' m ’ a * e w^sin 6 dw d0 dcp 


n£Tl nl T ^ 

Q(v) = \ \ \ vw^cos Be 

Jq Jq J o 


= 2tt 




o o 

w cosh^2j3j2 v w) w sinh^^gVw) 


%2 


2 %2 v 


”^12 w 
e 1Z dw 


, J» f "Uosh(2r C ) - cW£«le-C 2 d C 
/3 12 3J 0 L 2r J 


(B5) 


where, by definition, £ = ft 2 */ 2 w and r = ft^^v. ^ follows from equation (A15) 
that equation (B5) can be written as 

Q(v) = - ^ f ? 3 e" ?2 cosh(2r?)d? 

2 %2 ft 2 3 J 0 

= (r 2 - i)MM+ -Jl- f £e“ ?2 cosh(2r£)d? 

\ 2/ ft 2 ft 9 3j 0 


Mix) + 

2ir 

h2 

^12 3 

M(v) , 

27 r 

^12 

%2 3 


= (r 2 - + 27r - 1+ r C e”^ sinh(2r£)d£ 

V • ^12 0i, 3 L 2 J o ^ll 


(B6) 


where, in two instances, the factor £ exp (~? 2 ) was singled out for integration by parts. 
The remaining integral in equation (B6) can be transformed as follows: 

f e" ?2 sinh(2r?)d? = Je r2 T [e"^" r ) 2 - e"^ +r ^JdC 
Jq 2 Jq 


2 mL e ‘ s2<ls ■ C e ' s2<is ) 


e^fV^s 

Jn 


(B7) 


Equations (B6) and (B7), together with the expression for M(v) in equation (A18), yield 
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Q(v) = 


^ 12 ' 


r 2 2 
e r 




2 r 2 )*^0 



3 12 


r 2 2 
e r 


v- r» 1 1 „2n 

y (-D n+1 

n=o n'.(4n 2 - l)(2n + 3) 


(B8) 


Equations (Bl), (B3), and the first of equations (B8) provide the following equation 
for the drag or tractive force: 


D i2 " vrn 12 n l n 2°12 




2r 2 )'^0 



(B9) 


Thus, equation (B9) gives the drag 
or tractive force, per unit volume, 
exerted by the first gas on the 
second under the assumptions 
that each gas has a Maxwelliaq 
velocity distribution and that the 
molecules behave as hard spheres 
with the mutual collision diam- 
eter CTjgJ the symbol v denotes 
the bulk velocity of the first gas 
with respect to the second. 

The right-hand side of 
equation (B9) consists of a 
momentum vector vm^ multi- 
plied by the quantity 
defined by 



.1 1.0 10 100 


Dimensionless speed parameter, r 

Figure 6 .- Dimensionless ratio of Rjg to collision frequency N12 plotted as 
a function of the dimensionless speed parameter r. 


R 


12 - n l n 2 CT 12 


2 




+ 



2^0 



(BIO) 


which is plotted, as a ratio to the collision frequency given by equation (A19), in 

figure 6. Figure 6 shows that for r » 1, the quantity R^ is practically identical to 
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the collision frequency N^g, hut that, for r « 1, R^g is approximately^ 4/3 of N^g. 
These results can be verified analytically by observing that the expression in brackets in 
equation (BIO) is asymptotically equivalent to 


!Ki 


( r — o) 


(B11) 


and to 


1/2 

tt ' r 


(r - °°) 


(B12) 


As in appendix A, one is led to define a parameter k by 


|(x + | «-2) 7 = (i + ^)e- r2 + (2r ♦ 2 - - /j) 


r -s 2 
e ds 


(B13) 


or 


k = 


5/9 
2r 2 l 64 


(‘ + i?) 6 ’ r + ( 2r + r ' p)I 0 


r -s 2 
e s ds 


(B14) 


which is shown graphically in figure 2. For r « 1, k = 1, and for r » 1, k = 

X 2 o 

Equation (BIO) and the definition (B13) yield 


R 12 - 


n n rr 2 M 8 U . 2 „ r 2\ 


= ^(2u) 1/2 n 1 n 2 'J 


+fv 2 

W i2 5 


r 

, 1/2 


(B15) 


The corresponding expression for the drag or tractive force is 


D 12 " vm 12 R 12 


-|(2>) 1/2n ' 1 2Wl2 2 (^ + fv 2 ) 


/ 1 , «.,2\ 1/2 =: 


(B16) 


which is proportional to v for r « 1 and to v^ for r » 1. 

2 This approximation (applicable for r « 1) was first derived by Stefan (ref. 9). 
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RATES OF CHANGE OF MOLECULAR TRANSLATIONAL ENERGIES 
RESULTING FROM ELASTIC COLLISIONS BETWEEN MOLECULES 
OF TWO GASES WITH DIFFERENT TEMPERATURES 
AND DIFFERENT BULK VELOCITIES 

The time rate of change of the molecular translational energy per unit volume of 
the second gas, viewed from the second inertial frame, is derived in this appendix. 

Let a be a unit vector drawn from the center of molecule 1 of the first gas toward the 
center of molecule 2 of the second gas at the instant of collision. Let 0 W (0 ^ 9 W = 7t/ 2) 
be the angle which a makes with the relative velocity vector w = v^ - Vg immediately 
before collision. Then, 


— ► w ^ 

a = — cos 0 + A 

w w 


(Cl) 


where 


"A 


X • w = 0 


A = | A | = sin 6 = 

W (j 


12 


(C2) 


In the collision, molecule 2 receives the impulse 


I = 2m • o <7 = 2m cos 6^0 = 2m ^ cos 0 w ^w cos + w ^) (C3) 


The velocity Vg” of molecule 2 after collision is given by 


m 2 v 2" = m 2 v 2 + 1 


(C4) 


whence 


V ’ 2 - v 2 2 = —v 2 ' * + (— 

2 2 m 2 2 l mr 


(C5) 
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With the use of equations (C2), (C3), and (C5), the change in kinetic energy of 
molecule 2, resulting from the collision with molecule 1, can be written as 


1/2 

(C6) 

Since X • w = 0, by equation (C2), the decomposition of Vg into vectors perpendicular 
to w ^2 _l) an d parallel to w yields 



X = v 0 , cos (p. 


21 a 


12 


w 


(C7) 


where <p w is the angle Vgj^, X, that is, the angle between two half -planes emanating 
from the line of w which are parallel to the vectors Vg and a. 

The variables of integration are the six velocity components v^ . . . (rela- 

tive to the second inertial frame), the impact parameter b (or the corresponding 
angle 6-w)> and the angle cp w . The last of these enters the right-hand side of equa- 
tion (C6) only in the last term and therein, according to equation (C7), simply as 
cos < p w ; the integration of this term from <p w = 0 to 2 tt yields zero. Consequently, 
the time rate of change of the molecular translational energy per unit volume of the sec- 
ond gas, viewed from the second inertial frame, is simply 



where, by definition, 


S(v) 



- m 12 2\ “% w r¥2 ^ , 

w + w |we dv. dv„ 

m 0 I 1 2 


(CIO) 


32 


APPENDIX C 


With the change from the variables VpV 2 to the variables V,w defined by equa- 
tions (A4), equation (CIO) becomes 



w 


I we 


■(01+/3 2 ) v2 -0i2(w-vy 


dV dw 


(Cll) 


Since the term of the integrand of equation (Cll), of which V • w is a factor, can 
be resolved into three odd functions of the three components of V, its contribution to 
S(v) is nihil. Hence, upon recalling the first of equations (B3) and carrying out, for the 
other part of S(v), the integrations with respect to the components of V, equation (Cll) 
can be written 


S(v) = 


^1 + 02 


3/2 -0 
e 




— \w(v) + — P(v) 
02 / 02 


(C12) 


where, by definition 


W(v) = £ 


w^e 


-/3j 2 (w^-2v • w) 


dw 


(C13) 


Changing from the rectangular Cartesian coordinates w ,w ,w to spherical coord- 

x y z 

mates w, 0,<p, where the half-line whose co-latitude is 0 = 0 has the direction of v, 
and integrating over the intervals 0 S (p ^ 2n and 0 S 0 S n yield 

W(v) = f ? 4 e' ?2 sinh(2rC)d? (C14) 

/3 12 3 r J 0 

where, by definition, C = j3^ 2 and r = ®y singling out the factor 

C exp (-? 2 ) and integrating once by parts, a comparison of the results with the third of 
equations (A15) and the first of equations (B6) shows that 


W(v) = -2- M(v) + Q(v) (C15) 

012 

The second of equations (A13) and of equations (B3) together with equations (C12) 
and (Cl 5) give 


i 
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Furthermore, equations 


S(v) = - ^N(v) + ^ P(v) 


^12 


m< 


hj 


Ulr 


(Al), (A3), (the second of) (Bl), and (C16) yield 


2m 12 /m 12 <S 12 \ 

hf 12 


m 


12 


m 2 


vD 


12 


(C16) 


(C17) 
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DEFINITION OF A BINARY TEMPERATURE FOR 
CHEMICAL -KINETIC CONSIDERATIONS 

For considerations of the kinetics of chemical reactions, it is of interest to define 
a binary temperature T ^ which is a measure of the amount of translational energy 
available for chemical reaction in bimolecular collisions. For bimolecular reactions 
within a species, that is, when the reacting molecules are chemically identical, the tem- 
perature of the species serves a similar purpose and no newly defined special tempera- 
ture is required. 

For definiteness, let all velocities be measured, unless otherwise indicated, with 
respect to the second inertial frame. If, immediately before collision, molecule 1 has 
velocity w = v^ - Vg relative to molecule 2, the relative molecular speed along the line 
of centers is w * o = w cos 0___. (See appendix C.) This speed is the sum of the speeds 
w^ and w^ with which the two molecules approach their mutual center of mass; 
thus, 


w al + w ct 2 = w cos °w 


m l w crl = m 2 w a2 


The corresponding total kinetic energy is 


E a = \ m l w al + | m 2 w a2 = \ = | m 12 w2 f 1 “ ~ 2 

1 a 12 


and its mean value per collision is 


(Dl) 


(D2) 


1 (h^2) 3/2 C -% W 1 2 -^2 V 2 2 1 2C° 12 ( b 2 \ 


2irb db dv^ dVg 


(D3) 


1 2 

The inner integral (with respect to b) is equal to ^ ?rcr i2 ’ With the aid of equations (Al) 


and (A3), equation (D3) can therefore be written 

1 m 12 C 3„"^l w l t, " /3 2 v 2* 


E„ = 


a N(v) 4 


■f 


w e 


dv l dv 2 


(D4) 
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Transforming to the variables V, w by equations (A4) and carrying out the integrations 
with respect to the three components of V yield 


- _ m 12 / i t \ 3 / 2 ~Pl2 v2 W(v) 

a 4 ( /3 l + ^2/ N(v) 


(D5) 


where W(v) is given by equation (C13). Equations (A13) and (C15) can be used to write 
equation (D5) in the respective forms 


E 

a 4 M(v) 


m 


12 


2 , Q(v) 

P 12 M(v) 


(D6) 


Equation (D6), together with equations (A18), (B6), and (B7), yields 


E = 

" 4 * 12 


— + r 2 
2 


ir 

r J 0 


r -s 2 
e s ds 


-r" 


+ (l + 2r 2 )i^ e" s ds 


whence, for r « 1, 


(D7) 


E =l^fl + 2 r 2_Ar 4 + ...\ 
CT 2 p 12 \ 3 45 J 


(D8) 


and, for r » 1, the asymptotic formula 


E 


1 m 12/r 2 


2 


% 




(D9) 


is applicable. The content of equations (D7), (D8), and (D9) is conveniently expressed by 


° 2 /3 12 V 3 ) 


(DIO) 
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where 





(Dll) 


The parameter co, which is given by equation (Dll) and is plotted as a function of the 
dimensionless speed parameter r in figure 3, has the asymptotic values 1 and 3/4, 
respectively, as r — 0 and r — °°. 

It follows from equation (D8) that, for v = 0 and = Tg = T, 

I 0 = kT (D12) 

This result and either equation (D7) or preferably its equivalent, equation (DIO), provide 
a basis for defining a binary temperature T 12 which is a measure of the amount of 
translational energy available for chemical reaction in bimolecular collisions, that is, 


kT lH^H^ 2 ) <“M) 

or, to afford a better comparison with equation (A28), 

n kT i o o kT -i n kT « , « 

1 l+± 2 + co y 2 (D14) 

2 m^2 2 2 n^ 2 

When the comparison is made with equation (A28), it should be recalled that co and a 
are roughly equal to unity. Hence, T^ is approximately equal to another binary tem- 
perature T^ defined by equation (7) on the basis of the form of equation (A28); the two 
binary temperatures are, in fact, asymptotically equivalent, as r — 0. Similarly, they 
are asymptotically equivalent (as r — 0) to the expression for T^ given by equa- 
tion (D14) with w taken as unity. 

In summary, the preceding analysis suggests that the binary temperature T^ 
defined by equation (D13) may be a significant and useful parameter in theoretical con- 
siderations of reacting gas mixtures which are not in dynamical equilibrium. 
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NUMBER OF BIMOLECULAR COLLISIONS RESULTING 
IN CHEMICAL REACTION 


Assume that a fraction P of the bimolecular collisions, involving one molecule of 
each gas, for which the line-of-centers velocity component w • a (immediately before 
collision; see appendix C) lies in a small interval of positive real values, leads to chem- 
ical reaction. For simplicity, it is assumed that P is a step function, that is, 


P = 



= Constant 


(w • a < w) 
(w • a > w) 


(El) 


where W is the activation speed; the latter is related to the activation energy E a 
(see eq. (D2)) by 


E a =| m 12 * 2 


(E2) 


The quantity P 12 is the probability factor or steric factor for the reaction. 

In this appendix, the number of bimolecular collisions, involving one molecule of 
each gas and leading to chemical reaction, is calculated. This number will be denoted 
by N^g'- The sevenfold integration involves the impact parameter b and the 
molecular velocity components, that is, or, by appendix A, ^b,V,w) or 

(b,V,w, 6, cp). The last set will be employed here to express N^ 2 ' . 

It follows from equation (El) that the integration with respect to w need extend 

only from w = W to w = 00 ; for, when w < W, it is certain that w • a = w cos < W 

w 

and, therefore, P = 0. For given w > W, the condition P = P^ 2 is, by equation (El), 
equivalent 
tion (C2), 

= a 12 (l-cos 2 0 w ) (E3) 


equivalent to w • cT > W or w = w • a > W or IS cos # w > W/w. Since, by equa 


b = a l2 sin 0 w 


1/2 


1 1/2 


it follows that P = P^ 2 if and only if 
Let the quantity to the right of the inequality be denoted by b^(w). 


(E4) 
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Consequently, the techniques of appendix A and reference 13 can be used to write 


N 12' ° P 12"l"2 — i 7 fl t jV 


(E5) 


where 


n2v p Tt "(^i + / 3 2) v 2 " /3 12(^ _ ^) 2 f b i(w) 9 

= Jq Jq we 2irb dbw 2 sin 9 dw d0 d cp 


(E6) 


and 


(w - v ) 2 = w 2 - 2vw cos 9 + v 2 

Integrating with respect to b, ip, and the three components of V yields 

,2 


(E7) 


N ' — 2ir^/^P nn/3 2 / 2 2 0 ^12 V T 

w 12 _ Zir p 12 n l n 2*12 a 12 e h 


(E8) 


where 


T r°° / O 2^ -^io(w 2 -2vw cos e) 

I 2 = \ \ w(w 2 -W z je 12 sin 6 dw 6.9 


(E9) 


Now 
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2 - w 2 ~^ 12 w2 . 


W fi 12 v 


sinh^2j3 12 vw^dw = — ~ J (? 2 - z 2 )e ? sinh(2r?)d? (E10) 


where, by definition, r = z = /3 12 1//2 W, and $ = /3 12 i/z w. Singling out the 

factor £ exp (-? 2 ) for successive integrations by parts yields 

2 2 

J ? 2 e“^ sinh(2r£)d£ = - |-e“^ Jj sinh(2r£) + r cosh(2r£)] 
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e ^ sinh(2r£)dt; + Constant (Ell) 
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With this result, equation (E10) becomes 




sinh(2rz) + r cosh(2rzj] 


+ 






sinh(2r£)dO 


(E12) 


The changes of variables employed in equations (B7) can be used to reduce the integral 
in equation (E12); thus 


I 2 = 


-z 


2/3 


12 


[z sinh(2rz) + r cosh(2rz)j + e r ^] 


r 2 + - 
2 



(E13) 


Equations (E8) and (E13) finally yield the following expression for the number of 
bimolecular collisions, involving one molecule of each gas and leading to chemical 
reaction, 



12 n l n 2 cr 12 





(‘ 


+ 




(E14) 


which can be evaluated, for various combinations of z and r, with the aid of error- 
function tables and an asymptotic expansion of the complementary error function. 

For r = 0, equation (E14) reduces to 


N 12' * P 12 


2n l n 2°12 


'(i) 


1/2 


-z ‘ 


(E15) 


which is the product of the steric factor, the collision frequency (see eq. (A21)), and the 
factor exp(-z 2 ). The last factor can, by equation (DIO) with r=0 and equation (E 2), 
be written 
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(E16) 


40 



APPENDIX E 


This equation reduces, for = Tg = T, by equation (D12), to 

e - z 2 . e ' E a/ kT 


(E 17) 


which is the well-known Arrhenius factor of reaction-kinetics theory. On the other hand, 
if r * 0 and z = 0, equation (E14) reduces to = p i2 N 12 ( cf * e Q* (A21)). 

The observations of the preceding paragraph suggest the following approximate 
form for the more general but unwieldy equation (E14): 


N 


12 


= P 12 N 12 e 



(E18) 


where N^g is the collision frequency given by equation (A21) and T^g is the binary 
temperature given by equation (D13). In terms of the dimensionless parameters r 
and z, the approximation (E18) appears as 
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12 


= P 


12 2n l n 2 a 12 



+ 


2 

3 



(E19) 


The ratio of this approximate value of N^g' to its exact value (eq. (E 14)) is given in 
table I for several combinations of r and z. 

Since the approximation (E19) approaches zero too slowly as z approaches °°, 
the right-hand sides of equations (E14) and (E19) have been equated and solved for to. 
This value of to, denoted by to', which makes equation (E19) exact, is plotted as a 
function of r for various values of z in figure 4. 

A comparison of figures 3 and 4, in conjunction with the form of equation (E19), 
shows that, for z S 1, equation (E19) with to = 1 provides a better approximation for 
Ni 2 * than does equation (E19) with to given by equation (11) (plotted in fig. 3). 
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